Phonon-limited transport coefficients in extrinsic graphene 
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The effect of electron-phonon scattering processes over the thermoelectric properties of extrinsic 
graphene was studied. Electrical and thermal resistivity, as well as the thermopower, were calculated 
within the Bloch theory approximations. Analytical expressions for the different transport coeffi- 
cients were obtained from a variational solution of the Boltzmann equation. The phonon-limited 
electrical resistivity p e - P h shows a linear dependence at high temperatures, and follows p e - P h ~ T 4 
at low temperatures, in agreement with experiments and theory previously reported in the litera- 
ture. The phonon-limited thermal resistivity at low temperatures exhibits a ~ T dependence, and 
achieves a nearly constant value at high temperatures. The predicted Seebeck coefficient at very 
low temperatures is Q(T) ~ — ir 2 k%T/(3 e Ef), which shows a n -1 ^ 2 dependence with the density 
of carriers, in agreement with experimental evidence. Our results suggest that thermoelectric prop- 
erties can be controlled by adjusting the Bloch-Griineisen temperature through its dependence on 
the extrinsic carrier density in graphene. 
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PACS numbers: 65.80. Ck, 72.80.Vp, 63.22.Rc 

Recently, Efetov and KimA reported experimental mea- 
surements of the phonon-limited electrical resistivity of 
graphene samples, by achieving extremely high carrier 
densities by means of an electrolytic gate, in order to 
minimize the effect of other scattering mechanisms and 
to be away from the neutrality (Dirac) point. They ob- 
served that, at low temperatures, the electrical resistivity 
displays a ~ T 4 dependence, whereas at higher temper- 
atures it follows a linear ~ T trend, in agreement with 
an earlier theory by Hwang and Das Sarma£. The later 
based their analysis on the Bloch theory^, which ne- 
glects phonon drag effects and umklapp processes, and 
modeled the interaction between electrons and longitudi- 
nal acoustic phonons by the deformation potential ap- 
proximation. Moreover, they calculated the electrical 
resistivity* within the relaxation time approximation, 
which is valid for clastic scattering processes, a condi- 
tion not strictly satisfied in electron-phonon scattering. 

In this paper, phonon-limited transport coefficients in 
graphene are studied, by considering the presence of tem- 
perature and voltage gradients. Under similar assump- 
tions as in Ref. (— ), that is, Bloch theory and deforma- 
tion potential approximation for the interaction between 
electrons and longitudinal- acoustic phonons, analytical 
expressions for the transport coefficients are obtained 
from a variational solution of the Boltzmann equation^. 
This method, at the lowest-order approximation, pro- 
vides equivalent results to the relaxation time approxi- 
mation and, at higher orders, converges to solutions of 
the Boltzmann equation which are exact to an arbitrary 
precision, within the limitations of the set of variational 
functions chosen. As has been pointed out before*, the 
semiclassical approach based on the Boltzmann equation 
cannot capture the subtleties of electronic transport in 
graphene right at the neutrality point. This limitation 
can, at least in part, be attributed to the inhomogene- 
ity in the charge distribution experimentally observed in 
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FIG. 1: (Color online) Electrical resistivity calculated from 
Eq. (|31|) (lines), compared with experimental data (symbols) 
reported in Ref.(— ) at different carrier densities n in units of 
W 17 m~ 2 . Inset: Electrical resistivity at zero temperature, 
from data in Ref. (— ) (circles), at different carrier densities n. 
Also shown is the linear fit (dashed line) prr=o) = 14.276 + 
266.17 x 10 17 n _1 (fi) explained in the main text. 



graphene samples at the neutrality point 6 . However, at 
high carrier densities the Fermi level is well above the 
Dirac point, and the Boltzmann equation then provides 
a correct description. 



I. ELECTRON-PHONON INTERACTION 

As is well known from tight-binding calculations, elec- 
tronic transport in graphene displays relativistic fea- 
tures. This is because pristine graphene, free of impu- 
rities and defects, behaves as a metal only in the so- 
called Dirac points where the conduction and valence 
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bands touch. There exists two non-equivalent such points 
within the first Brillouin zone, defined by the vectors 
Kx,2 = ±— [l/\/3, l], with a the direct lattice param- 
eter. In the vicinity of each of these points, or val- 
leys, it is shown that the dispersion relation is linear, 
E n (8'k) = ±ftvi?|<5k|, for <5k = k K., ; . In each val- 
ley K l; , the envelope electronic functions related to this 
dispersion relation can thus be obtained as the pseudo- 
spinor solutions of an effective two-dimensional Dirac's 
equation, 

H^A ±] = v F a ■ pi{ ±] = Sf^f ) • (1) 

Here a = (a x ,a y ) are the Pauli matrices, we have re- 
defined k = 5k for each valley, and the energy eigenvalues 
for positive (negative) helicity are E ( k ±} = ±hv F k. The 
corresponding eigenvectors, normalized by the total area 
A of the sample, are given by the pseudo-spinors 

with tfk = arctan(fc y /fc x ). 

In the language of second quantization, the conduc- 
tion electrons field for k in the vicinity of the K^2 = 
±^ [l/\/3, 1] valley is 

^w = E^ +) ( r )^ (3) 

k,tT 

with Cko-r; Fermionic operators [cko-ij, cj c , (T , ?) ,]+ = 
^kjk'^o-jCT'^r),?)' , that destroy (create) a conduction elec- 
tron with momentum Kk in the vicinity of the K, ( val- 
ley, with spin component a = {t, 4-}- We model the 
electron-phonon interaction by the deformation poten- 
tial approximation 2 -, which is a reasonable one when 
the Fermi surface possesses spherical symmetry^. The 
two-dimensional version of this is a Fermi "circle" in 
graphene, which corresponds to the intersection of the 
Dirac cone and the constant energy plane defined by the 
Fermi energy at eu = Ep. In the deformation potential 
approximation, the operator representing lattice dilation 
in the second quantization language is 

A(r) = iJ2 (2A Pm u Jq y 1/2 (q • e L ) (a\e^ r - a q e-^')(A) 
q 

Here, to q = v s q is the frequency for longitudinal acoustic 
phonons in graphene, a q is the creation operator for longi- 
tudinal acoustic phonon modes, and is the correspond- 
ing polarization vector for these modes. The electron- 
phonon interaction Hamiltonian, in the deformation po- 
tential approximation, is obtained from 

H?* ph = * [d 2 rJ2 $(r)DA(r)^(r) (5) 

tj=l,2 



Notice the presence of the factor cos(#/2) = cos([(fk — 
(pk']/2) which arises from the inner product of the two 
pscudo-spinor functions u£, • defined in Eq. ([2]). This 
factor, which would be absent from the inner product 
of non-relativistic scalar Schrodinger wave-functions, is 
a fingerprint of the relativistic features of electrons in 
graphene. The more obvious consequence of its presence 
the suppression of backscattcring (8 = 7r), a phenomenon 
observed in graphene and directly related to so called 
Klein tunneling^— in the context of relativistic quantum 
mechanics. 

We shall neglect umklapp processes, and thus we have 
G = for normal processes. In this case, scattering 
events satisfy energy and "crystal momentum" conserva- 
tion, ey = tk^fajJq and k' = k±q. From the latter equa- 
tion, one obtains the condition q 2 = k' 2 + k 2 — 2kk'cos(6), 
with 8 the scattering angle. Since the relevant scatter- 
ing events occur within a thin layer near the Fermi sur- 
face, then k ~ k' ~ kp, which then implies \sin{8/2)\ = 
q/2k F . In a Debye model approximation, this condi- 
tion restricts the range of lattice momenta involved in 
electron-phonon scattering, q < min{2/cp, to}, with qp> 
the radius of the Debye circle (for a two-dimensional sys- 
tem). In geometrical terms, the most restrictive condi- 
tion depends on whether the radius of the Fermi sur- 
face (a circle centered on each Dirac point for graphene) 
is larger or smaller than the Debye circle. When writ- 
ten in terms of the frequency of the longitudinal-acoustic 
phonon modes involved u> q = v s q, the former condition 
reads Huq/kg < minjOsGj @-d}j with 0£> the Debye 
temperature, and &bg = 2(v s /vp)Ep/kB the Bloch- 
Griincisen temperature. In the experimental setup pre- 
sented in Ref.(— ), the carrier densities involved are such 
that the Fermi surface is contained within the Debye 
circle, and therefore the Bloch-Griineisen rather than 
the Debye temperature imposes the characteristic en- 
ergy scale for phonon scatterers. Notice that the analy- 
sis of the scat tering angle ab ove also yields the relation 
cos(0/2) = v 71 - (<7/( 2 M) 2 - Therefore, after Eq. © 
the matrix element for the electron-phonon interaction 

H l e - ph = % Ek,q,„,<7 M q (oq ~ «-q) 4 +q ,^ 6 k CT7) IB 

1 /2 

M q = -l (h/2p m tu q ) 1/2 Dq (l - (q/2kpf^ . (6) 

Here, D is the deformation potentiali^ coupling con- 
stant for graphene, and p m = 7.6 x 10 -7 kg/m 2 the sur- 
face mass density. 

II. BOLTZMANN EQUATION 

We introduce, in the usual form, the deviation of the 
distribution function from equilibrium by the expression 
fk = — Xkdfk/dtk- Let us consider the electron- 
phonon interaction as the sole scattering mechanism. 
Therefore, after Eqs.(5|H]) we should consider processes 
of the form k + q <=± k'. Under Bloch's approxima- 
tion, which is to assume that the phonon system is in 
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quasi-cquilibrium, the linearized form of the Boltzmann 
equation when both an external electric field and a tem- 
perature gradient are imposed can be expressed as 



-v k -VT^I-v k 



dT 
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k',q 



Here, g n = g a = 2 are the spin and valley degeneracies, 
respectively, while the transition rate for electron-phonon 
scattering processes is 
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x</°(l-A ,) 



q |M,|^(e fc - e k , + huj q ) 



(8) 



Here f% is the Fermi-Dirac distribution for electron states 
with momentum k, and n° the Bose distribution for 
phonon states with momentum q. All vectors are two- 
dimensional, with G a vector of the reciprocal lattice. 

For extrinsic graphene, the Fermi level Ef = hkFVF is 
located above the Dirac point, and depends on the den- 
sity of carriers n by the relation kp = \/irn. In Rcf . (— ) , an 
experimental method is presented which allows to control 
the carrier density, and hence the Bloch-Griincisen tem- 
perature. As discussed in Refs.(— this in turn becomes 
a way to control the phonon-limited electrical resistivity. 
In this work, we will show that also the electronic com- 
ponent of the thermal conductivity, as well as the ther- 
mopower (Seebeck coefficient) which determines thermo- 
electric effects, can be controlled in a similar way. We 
shall use the variational method to calculate the phonon- 
limited transport coefficients which, in contrast with the 
relaxation time approximation, it has the advantage of 
avoiding the assumption of quasi-clasticity in the scat- 
tering processes. 



III. COUPLED CHARGE AND HEAT 
TRANSPORT 

For a system with thermal and electrostatic poten- 
tial gradients, the macroscopic entropy production rate 
is given in terms of the heat U and charge J fluxes by^A 



d t S = J-E/T + U- V(l/T), 



(9) 



The Onsager's theorem^ indicates that the fluxes and 
generalized potentials are then linearly coupled, 
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(10) 



The thcrmopower Q can be expressed in terms of the 
coefficients L a p by considering the case J = in Eq. 
(fTU|> . to obtain E 
that Q = —LetL 



J EE- 



LetL ee VT . Hence, one concludes 
The thermal conductivity is ob- 



tained from the second equation, U = —{LteL ee Let~ 
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FIG. 2: (Color online) Phonon-limited thermal resistivity 
[K e - P h]~ , calculated from Eq. (|34[) at different carrier den- 
sities n, in units of 10 17 m -2 . Shown in the inset is the cor- 
responding thermal conductivity, including the contribution 
from Coulomb elastic scattering with impurities, and normal- 
ized with a nominal packing thickness of 3.4 A. 



Ltt)VT, which implies n = LteL ee L E t — Ltt- If. 
on the other hand, one considers the case VT = 0, then 
the relation J = L EE ~E is obtained. Thus, it is con- 
cluded that the electrical resistivity is p = L EE . The 
variational method^ is based on the principle of maxi- 
mal entropy production^, and provides an iterative but 
virtually exact procedure to solve the Boltzmann equa- 
tion, and therefore to calculate the transport coefficients, 
within the limitations of the set of variational functions 
chosen. The principle is based on equating the macro- 
scopic entropy production rate, as expressed in terms 
of the macroscopic currents Eg . (fit)]) . with the entropy 
production rate due to microscopic scattering events^, 
dtS = S SC att, as defined from the linearized Boltzmann 
equation scattering terms, 



knT 2 



J2 {Xk-Xk'} 2 ^] 



k' 
kq 



(11) 



k,q,k' 



The electric current J and the heat current U are given 
by the expressions 



9v9<7^2i~e)-v kX k 



M 



u = 9ri 9a v ^ - ^)^ k fr (12) 



The technique is based on expanding the distribution 
function deviation from equilibrium Xk = Y)j otj&iO*) 
in terms of a set of variational functions of the form 
4>i(k) = (efc — EpY^iu ■ k). For a minimal set of trial 
functions 0i(k) = k ■ u and 02 (k) = k • u(ek — E E ), we 
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obtain the charge and heat currents 



J* =(-e)ft l .9aVv k (j!) I (k)^- 
, Cefe 



Eq.Ql 



df° 



Ui = ff^ffo- ^ Vk(efc -E F )(f>j (k) -g— (13) 



Here, g, ; = g a = 2 are the valley and spin degeneracies 
in graphene, respectively. The coefficients Pij = Pji are 
defined as follows 



ij 

Let = -T-^J^P- 1 ]^^- 

y 

Ltt = -T~ 1 Y^'Ui[P~ 1 ]i3^j 



(19) 



13 k B T 



P,, = 



£ {<Mk) - ^(k')}{^(k) - ^(k')}^ q 



k.k'.q 



(14) 



In terms of the variational solution, after Eq.(|T2j) the 
electric and heat currents arc written as 



Therefore, we proceed to calculate the currents and 
Pij coefficients. Let us first consider the current Ji. 
The group velocity in Eq. (fT3f is given by Vk = kvp , 
with k = k/fc = (cos(y), sin(y)). The discrete sums 
are treated in the usual quasi-continuum limit as two- 
dimensional integrals, ^ k — » A(2ir)~ 2 J cPk, and the 
integral is calculated by using the change of variables 
d 2 k = dtpede / (hv f) 2 , 



J = ^ ajJj 

i 

U = ^a ?; U 4 



Ji 



(15) 



„Ae f E F \ 2 
nh \hvp J 



(20) 



whereas the macroscopic d t S and microscopic S scat t en- 
tropy generation rates in Eqs. (|9lll[) are expressed by 



d t s = J2 a i ■ E - T_2u * ' VT ] - 



For the current Ui, we consider the quasi-continuum 
limit for the sum over wave- vectors in Eq. ([13j) . with direc- 
tions k = k/fc = (cos(ip) , sm(tp)) . Using the integration 
variables d 2 k = dipdee/ (hvp) 2 , we have 



Sscatt = T 1 Qua j Pij. 



(16) 



A 2-kA (k B T) 2 E F 



The variational problem is to find the set of coefficients 
{ay} that maximize Sscatt, subject to the constraint 
dtS = S sca tt, which is enforced by a Lagrange multiplier 
A, 



3 h (hv F ) 2 

Similarly, for the current U2 we obtain 
Av F 



(21) 



U 2 



ir 2 (hv F ) 3 7_ 7r 



difkk-u I dee 2 {e-E F f 1 
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P 



J.-E 

T 



VT 

rp2 



^A (k B T) 2 / E F 

3 h \Hvf 



(22) 



(17) 



The solution to this constrained variational problem is 
given by A = 1/2, and 



(18) 



Here, [-P -1 ] is the inverse of the matrix [P]ij = Pij whose 
elements are the coefficients Pij defined by Eq. (fT4]) . Af- 
ter substituting this solution into the expressions for the 
electric and heat currents, we obtain explicit analytical 
formulas for the transport coefficient tensors defined in 



Finally, it is straightforward to verify by inspection that 
J 2 = (-e)Ui. 

Let us now consider the coefficient Pn. From Eq.©, 
we first notice that the transition probability rate 
imposes the condition q = k' k. Wc define the 
directions k = (cos(<//), sia(^')) and fc = (cos(<// — 
0),sin(y>' — 9)), with 8 the scattering angle. As before, 
in the quasi-continuum limit for the sums over wave- 
vectors, we choose the integration variables d 2 kd 2 k' = 
dip'd9de'e'dee/(hvp) 4 , to obtain 



Pu 



(ksT)-^ 2 



dee / de'e' 



{whvp) 4 

I del V(q-«) 2 ? k k q (23) 

J —7T J — 7T 
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In Appendix A it is shown that the angular integral 
/ dy'(q- u) 2 = 7r<7 2 , with q = 2k F sin(9/2). Hence, sub- 
stituting the explicit expression for the transition proba- 
bility rate, Eq. (j2"3")) reduces to 



Tr 2 h(hv F ) A 

x / dee(e+^)/ e °(l-/° ) (24) 



The energy integral is evaluated using the technique pre- 
sented in Appendix B. Using the change of variables 
dO = k F l [l - {q/{2k F )) 2 }- l ' 2 dq, we obtain 



(k B T)~ 1 A 2 k F HD 2 f 2k - f 4 
Pl1 = 2 71 — ^2 / "9 9 



,v F 



T 



e 



BG 



2„4 



7T q 



<£ 
2 



[1 



g 
2k F 



il/2 



X(l - e - h ^J( k BT)yl^ q /(k B T) _ 



(25) 



Eg. (1231) is written in terms of the functions J p (z) defined 
in Appendix C, 



Pn 



hD 2 A 2 E F (k B Tf 

TT 2 p m (flV s ) 5 (hv F ) 3 

T 



MT/Q BG ) + A 



v F 



2 (ir 2 1 
Wbg / I o 2 



(26) 



Let us now consider the coefficient P22, which after sim- 
ilar manipulations as before, can be expressed by the 
integral form 



1 a2 r* 



P? 



dO 



dip' / dee / de'e 



(k B T)~ A 
(irhv F ) 4 

x[(q -u)(e-E F ) + fiw,(k' • fi)] 2 P& 



(27) 



After expanding the square in Eq. (|27|) . it is shown in 
Appendix A that the angular integrals give, J dip'(q ■ 
u) 2 = nq 2 , JV(k' • u) = irk' 2 , and /<V(q • u)(k' • 
u) = (7r/2)[g 2 + fc' 2 — fc 2 ]. Then, after similar algebraic 
procedures as described previously, Eq. (|2"7|) reduces to 



hD 2 A 2 E%{k B T) i 

Tr 2 p m (hv s ) 3 (hv F y 



An 2 , 
1 + -3- (1 + 12 




T \ 2 
7 —) MT/@ BG )--{1-Q 



Bbg/ 3 \v F 

v £)MT/e BG ) + ^{i 



Js(T/e BG ) 



(28) 



At last, we consider the coefficient Pi 2. After our defini- 
tion Eq. tfbfj) . we have 



Pu 



dO 



dip' / dee / de'e 



(k B T)~ 1 A 2 

x [(q • u) 2 {e - E F ) + hu q (k' • fi)(q • u)]P£ q (29) 



As shown in Appendix A, the angular integrals are given 
by J dip'((\ ■ u) 2 = nq 2 , and J dip' (q ■ u)(k' • w) = 
(7r/2)[g 2 + A:' 2 — fc 2 ]. Therefore, using the procedure de- 
scribed previously, the integral in Eq. (|29)) is calculated 
and expressed in terms of the Jp(z) functions defined in 
Appendix C, 



P12 = 



HD 2 A 2 E 2 F (k B T) 4 
ir 2 p m (Hv s ) 3 (hv F ) 5 
8ir 2 ( T 



J A {T/Q BG )-^(-^- 
vf \®bg 



x J 5 (T/Q BG ) + 



Vs 
V F 

Vs 
V F 



MT/<d BG ) 
3 



6bg 
1 



MT/e 



BG) 



j2 MT/Q BG )\ + - (30) 



T 



e 



BG 



{J 7 (T/Q BG ) - 2tt 2 J 5 (T/Q bg )} 



IV. ELECTRICAL RESISTIVITY 

As discussed in Section III, the electrical resistivity 
is obtained by setting VT = in Eq. (fTQ|) . to obtain 
p = Lfig. By direct substitution of the coefficients and 
currents obtained in Section III, we obtain after Eqs. fTi?)) 
that the leading term which defines the electrical resis- 
tivity due to electron-phonon scattering processes is 



p e -ph{T) = pa 



"BG / 3 \V F 



T 



e 



BG 



MQbg/T) + —jMQbg/T) 
2it z 



(31) 



Here, po = 8 D 2 k F /(e 2 p m v s v F ) is a coefficient with 
dimensions of a two-dimensional resistivity (O), and 
Q BG = 2(v s /v F )E F /k B is the Bloch-Griineisen temper- 
ature. The functions Jp(@ BG /T) and their asymptotic 
properties are defined in Appendix C. Eq. (f3"Tj) depends 
on the dimensionless parameter T/Q BG , a feature pre- 
viously obtained^ within the relaxation time approxima- 
tion. The functions J P {Q BG /T), defined in Appendix 
C, have the property Jp(oo) = p\(,{p), with C,(p) the 
Riemann zeta function. Therefore, the low-temperature 
(T <C Obg) behavior of the electron-phonon contribu- 
tion to the electrical resistivity is 

Pe- P h(T) = po 4!C(4)(T/9 BG ) 4 + o(T/e BG ) 6 (32) 

At high temperature (T > Q BG ), from Eq. lpjT]) 
we find the asymptotic limit p e - p h{T) = {1/3 — 
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1/10... }p T/Q B G = (tt/8)p T/Q bg . This behavior is 
in agreement with experiments 1 - and earlier theoretical 
results based on the relaxation time approximation 2 -, thus 
supporting the application of the variational method in 
this case. In Fig.([T|) we compare the theoretical predic- 
tion with the experimental data reported in Ref. (— ). 
The experimental data are fitted to the expression 

p(T) = p (T =o) + Pe-ph{T), (33) 

with pe-ph (T) defined by our theory Eq. (|31|) , and /0(t=o) 
being the value of the resistivity at zero temperature. 
This parameter can be interpreted as the (temperature 
independent) contribution of impurity scattering to the 
total electrical resistivity. The impurity contribution 
to the electrical resistivity has been discussed 5 -^, and 
it is a sample-dependent property, since it is propor- 
tional to the surface concentration of impurities rii mp . 
It has been showr* 5 -^ that, for a short range scatter- 
ing potential V sca tt(^) = Vb<5(r), the in-plane resistivity 
is Pscatt oc Tii m pVQ, whereas for a long range Coulomb 
potentiate p sca tt oc 7ti mp rj _1 , this last case showing an 
inverse dependency on the carrier concentration n as well. 
Therefore, it is expected that the zero temperature resis- 
tivity exhibits a dependence on the carrier concentration 
of the form prx=o) = a + bn , with the value of the 
coefficients a and b depending on details of the sample, 
particularly on the concentration and distribution of im- 
purities and defects. We verify that this equation is in 
excellent agreement with the data reported in Ref. (— ), 
as shown in the inset of Fig. (TTJ). We extracted values 
a = 14.276 (fl) and b = 266.76 x 10 17 (flm 2 ), with linear 
regression coefficient r = 0.999. 

As seen in Fig. (TTJ), Eg. ([33)) provides an excellent fit 
to the experimental data. The temperature-dependent 
electron-phonon contribution to the electrical resistivity 
is fitted with just two free parameters: The energy pa- 
rameter D — 23.5±0.5 cV in the deformation poten- 
tial, and the speed of sound for longitudinal phonons 
in graphene v s = 24.0 ± 0.6 Km/s. It is remarkable 
that both values are in excellent agreement with inde- 
pendent estimations reported in the literature, particu- 
larly the speed of soun d 12 ' 13 . It is also relevant to notice 
that the ratio (v s /vf) 2 ~ 10 -4 is negligibly small, and 
hence terms proportional to this factor can in practice be 
neglected in the analytical expressions obtained for the 
electrical resistivity, as well as in other transport coeffi- 
cients discussed along this work. This also explains the 
quantitative success when applying the relaxation time 
approximation to electron-phonon scattering as reported 
in previous theoretical work 2 -, even though this scattering 
mechanism is not strictly quasi-elastic. 



V. THERMAL RESISTIVITY 

The thermal conductivity is obtained from the gen- 
eral expression Eq. (flU)) by setting J = 0, which as 
discussed in Section III leads to the expression n = 
LteL'^Let — Ltt- By direct substitution of the coeffi- 



cients and currents calculated in Section III, and neglect- 
ing terms proportional to (v s /vp) 2 ~ 10~ 4 , we obtain 
that the leading contribution to the thermal resistivity 
due to electron-phonon scattering is given by 



[He-ph(T)] 1 = 

x MQbg/T) - 



Po 

C T 
1 

' 2^2 



{T/® BG f + —{T/e BG y 



{T/Qbg) MQbg/T)\ (34) 



Here, Co = Tr 2 k B /(3 e 2 ) is the Lorenz number for the free 
electron gas. At very low temperatures (T -c ®bg), Eq. 
(|34|) predicts the asymptotic limit 



CoObg 4tt 2 



5 4! C (4)- 3 



e 



BG 



T 
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BG 



(35) 



It is remarkable that, according to Eq. (|35p. the phonon- 
limited thermal resistivity is linear at very low tempera- 
tures, in contrast with the typical ~ T 2 behavior theoret- 
ically predicted and experimentally observed in normal, 
three-dimensional metals 3 -"^. The thermal resistivity con- 
tribution due to elastic scattering with Coulomb impuri- 
ties has been calculated 5 -, and corresponds to [K SC att] _1 = 
3hT- 1 ul/(2n 2 k 2 B E 2 F ), with = n imp Ze 2 /(16e 2 eg). The 
total thermal resistivity can be estimated from the ex- 
pression [k(T)]" 1 ~ K-phiT)}- 1 + focttCT)]- 1 . In 
Fig.@ inset, we compare the value of the total three- 
dimensional thermal conductivity at room temperature, 
normalized by a nominal packing "thickness" of 3.4 A 
for the graphene layer. For numerical evaluation, we as- 
sumed Z = 1 for the impurities, and from the exper- 
imental values of the zero temperature electrical resis- 
tivity, we extracted an average impurity concentration 
n-imp = 1.3 x 10 15 m . For the experimental system 
reported in Ref.(— ), we estimated a dielectric constant 
e = 3.1 representing the average between the SiO? sub- 
strate and the PEO polymer electrolyte. As seen in the 
inset of Fig.©, the total normalized thermal conductiv- 
ity at room temperature is on the order of K ~ 400 Wm _1 
K _1 . We can compare this with the thermal conductiv- 
ity due to the phonon system, where the experimental 
valu o 14 ' 15 , in excellent agreement with a theory previ- 
ously reported by the present author—, is abou t 12 ' 14 ' 15 
K p h = 4 300 Wm^'K -1 at room temperature. These 
results suggest that, in most of the temperature range, 
the phonon contribution to the thermal conductivity in 
graphene dominates over the electronic contribution, as 
the present author has pointed out elsewhere^ 2 -. 



VI. THERMOPOWER 

As discussed in Section III, by setting J = in Eq. (fTU|) . 
we obtain that the thermopower (Seebeck coefficient) Q 
is given by Q = —Let/ Lee- By direct substitution of 
Eq. p^|) . we obtain that the leading contribution to the 
thermopower is given by the expression 
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Q(T) = - 



ir 2 k%T 




MObg/T) + % (^) J 5 (T/Q BG ) - 2 [^-^ J 6 {Q BG /T) 


3e E F 


1 + 


i^(_i_\ 2 

3 ^e BG y 


<Ji{® bg /T) - | (e^) 2 MT/& BG ) 



(36) 
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FIG. 3: (Color online) Thermopower at different electronic 
densities n, in units of 10 17 m -2 , calculated after Eq. (|36[) with 
v s = 24 Km/s as inferred from experimental data in Ref.(— ). 
Inset: Experimental data (symbols) for hole thermopower in 
Ref.(— ), compared with Eq. (|36p at the corresponding hole 
densities (no fitting parameters besides the value of v s — 24 
Km/s are used). 



VII. CONCLUSIONS 

We presented a semi-classical theory to calculate 
the phonon-limitcd transport coefficients in extrinsic 
graphene at high carrier densities. Our theory is based 
on a variational solution of the Boltzmann equation, 
and provides explicit analytical expressions for the var- 
ious transport coefficients, such as electrical resistiv- 
ity, thermal resistivity and thermopower. We showed 
that our analytical results are in excellent agreement 
with experimental data arising from two independent 
groups, particularly concerning the electrical resistivity 
and thermopower—. 

Our theoretical results suggest that, in principle, it is 
possible to control not only the electrical resistivity, but 
also the thermal resistivity and thermopower by control- 
ling the extrinsic carrier concentration in graphene. 
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Therefore, we find that the thermopower at very low 
temperatures (T <IC @bg) becomes 



Q(T) 



TT 2 k%T 



5!C(5) 



w F 4!C(4) V© 



T 



>BG 



o(T 3 X37) 



This result is in agreement with experimental data re- 
ported in the literature^. In particular, it has been 
observed^ that the thermopower in graphene depends 
on the carrier density as ti^ 1 / 2 , and shows a linear in T 
dependence at very low temperatures. Since the Fermi 
level is Ep = YmF^~wri, it is clear that Eq. ([36| obtained 
from our theory correctly reproduces this feature. The 
contribution to the total thermopower due to phonon- 
drag effects has been discussed in Ref.(— ), where it is 
shown that it displays a T 3 dependence at low temper- 
atures. Clearly then, the phonon drag is a negligible 
contribution to the total thermopower as compared to 
the diffusion component calculated in this work, which in 
agreement with the experimental data^ reproduces the 
correct linear in T dependence at very low temperatures. 
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Appendix A: Evaluation of angular integrals 

We consider the elementary integrals over angular ori- 
entations of the wave- vectors used in the main text. Let 
us define the unit vector u = (cos(a), sin(a)), with a an 
arbitrary (but fixed) direction along the thermal gradi- 
ent or the applied electric field. As in the main text, 
we define the scattering angle by k ■ k = cos(0), 
with k, k! the directions of the initial and final wave- 
vector after the electron-phonon scattering event takes 
place. Then, we write k — (cos(y'), sm(<p')) and k — 
(cos(< / 9 / + 9),sin(ip' + 8)), respectively. Using this nota- 
tion, we consider the following integrals: 



dip' [k ■ u 



r,l2 



+7T 



dtp' [cos (a) cos(<y5' + 6) + sin(a) sm((p' + 0)} 
7r[(cos(a)) 2 + (sin(a)) 2 ] = tt (Al) 
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f + TT 

dtp'[k' ■ ii} 2 

— TC 

dtp' [cos(a) cos(<^') 4- sin(a) sin(<^')] 

-7T 

= 7r[(cos(a)) 2 + (sin(a)) 2 ] = tt (A2) 

dip [k ■ it] [k ■ u] 

dip' [cos (a) cos(p' + 9) + sin(a) shx((p' + 6)] 

x [cos(a) cos(y') + sin(a) sin(<p')] 

dip' [(cos(a)) 2 cos(p' + 6) cos(p') 

+ (sin(a)) 2 sin((^' + 6) sin(y/) + sin(a) cos(a) 

x sin(2<// +0)]=n cos(0) (A3) 

Direct application of these basic identities leads to the 
following results 



+ 7T /* + 7T 

dp' [k ■ u] [k' ■ u] = kk! I dp' [k ■ u] [k' ■ u] 



irkk' cos(0) 



(A4) 



[ " dp'[ci-u} 2 = f "V[(k'-k)-u] 2 

= TT(k 2 + k' 2 - kk' COS(0)) = 7T(k' - k) 2 

(A5) 



irq 



dp' [q • u] [k' -u]= dp' [(k' - k) • u] [k' • u] 

-TT J —TT 



w(k' 2 - 2fcfc'cos(6>)) = tt — + 



q 2 k' 2 - k 2 



(A6) 



Appendix B: Evaluation of energy integrals at low 
temperature 

In several calculations throughout this work, after 
the change of variables x — (e — Ep) / '(fc^T) and z = 



hjj I '(fegT), we use the identity 



F(x) 



-dx = (l - 



< y + °° [G(x) - G(s - z)] (~^Pj dx (Bl) 



with G(x) = Jq F(x')dx' . This identity is straightfor- 
ward to prove using integration by parts, and noting that 
f°(x) - f°(x + z) = (1 - e- z )(l + e*)" 1 ^ + e~ x - z )- x . 

The right-hand side of the equation can be further 
evaluated using the Sommcrfcld expansion for the Fcrmi- 
Dirac distribution function at low temperatures, leading 
to the formula 



F(x) 



(1 + e*) (1 + e-*-*) 



dx 



F(x)dx+—[F (O)-F(-z)} 
6 



(B2) 



which is valid up to o{(ksT '/ 'Ep) 4 



Appendix C: Definition and properties of the 
functions J v {z) 



We have defined the functions 



2 lPy/1 - (x/z) 2 

(1 - e- x )(e x - 1) 



dx 



(CI) 



The asymptotic limit of these functions, for p > an 
integer, is given by J~ p (oo) = p\(,(jp), with the Rie- 
mann Zeta function. In particular, for p = 4 and 6 it has 
the values C(4) = tt 4 /90 and C(6) = tt 6 /945. 
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